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We report the surface electronic structure of niobium phosphide NbP single crystal on (001) surface by vac¬ 
uum ultraviolet angle-resolved photoemission spectroscopy. Combining with our first principle calculations, 
we identify the existence of the Fermi arcs originated from topological surface states. Furthermore, the sur¬ 
face states exhibit circular dichroism pattern, which may correlate with its non-trivial spin texture. Our results 
provide critical evidence for the existence of the Weyl Fermions in NbP, which lays the foundation for further 
investigations. 

PACS numbers: 71.55.Ak, 74.20.Pq, 74.25.Jb, 79.60.-i 


In the past few years, great progress has been witnessed 
in the study of quantum materials with non-trivial topological 
electronic structures. By introducing topological orders, insu¬ 
lators can be further classified into trivial one and non-trivial 
one, that is, topological insulators (TIs) dEl- Due to their 
unique physical properties and great potential in applications, 
TIs have been one of the central research subject of condensed 
matter physics in the last decade. 

The recent proposal that Weyl fermions can be realized in 
the so-called topological Weyl semi-metals (TWSMs) broad¬ 
ens the classification of topological phases of matter beyond 
TIslUCl. The low energy excitations in this new type of topo¬ 
logical quantum matter are described by the Weyl equation (SI . 
Thus, in the bulk, the conduction and valence bands disperse 
linearly cross pairs of discrete points (the Weyl points) along 
all three momentum directions. The Weyl points are associ¬ 
ated with a chiral charge that protects gapless surface states 
(SSs) on the boundary of a bulk sample. These topological 
SSs take the form of unclosed curves connecting the Weyl 
points of opposite chirality (H, leading to the existence of 
the unique Fermi arcs on the surface. Due to the novel phys¬ 
ical phenomena related to the Weyl fermions such as nega¬ 
tive magneto-resistivity, quantum anomalous Hall effect and 
non-local quantum oscillations, the search of TWSMs have 
attracted worldwide attentions lOUTlIl. 

In principle, by breaking either time reversal or spacial in¬ 
version symmetry of a Dirac semi-metal (DSM) to remove 
the spin degeneracy of the bands |[T5ti23]| . the otherwise de¬ 
generate bulk Dirac point will be split into a pair of Weyl 
points of opposite chirality, leading to the realization of the 
TWSM. Early predictions of the TWSMs focus on magnetic 
materials, e.g. R2lr207, HgCr2Se4, which naturally break the 
time reversal symmetry naia. However, due to the intricate 



FIG. 1: (Color online) The crystallographic and transport properties 
of NbP. (a) Illustration of body-centered-tetragonal structure, (b) The 
X-ray diffraction pattern of the (001) surface. The full width at half 
maximum of the (004) peak in the rocking scan is no larger than 
0.052°. (c) Temperature dependent resistivity curve at zero magnetic 
field, with RRR over 100. The inset presents the image of the as- 
grown single crystals, showing fiat and shiny surface, (d) Bulk and 
(001) surface Brillouin zone. Twelve pairs of Weyl points are de¬ 
noted by red and blue dots to represent opposite chirality. Note, the 
separation of each pair of Weyl points is exaggerated for clarity. 


magnetic domain structure of the samples, the experimental 
verification of these compounds still remains an open ques¬ 
tion. Alternatively, the non-centrosymmetric transition metal 
monoarsenides/phorsphides, TaAs, NbAs, TaP and NbP, have 
been proposed to be candidates of the TWSMs very recently 
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FIG. 2: (Color online) The measured and calculated surface states (SSs) of NbP on the (001) surface, (a) The photoemission intensity map 
integrated over [Ep - 10 meV, Ep + meV]. (b) The calculated FS which nicely reproduce the measured one. (c), (d) Three-dimensional 
intensity plot of the photoemission spectra around the X and F, respectively, (e), (f) The measured and calculated band dispersions along Y - 
M - X, respectively. The size of the red points in the calculated FS and band dispersions scale with the weight projected on the first unit cell of 
the P-terminated side, (g) Photon energy dependence of the momentum distribution curves integrated near Ep along high symmetry directions 
E - X and E - F, respectively. All the presented data [except those in panel (g)] were taken with circularly-polarized 120 eV photons at 12 K. 


GiEa. These compounds naturally break the inversion sym¬ 
metry, meanwhile avoid the complication of working with cor¬ 
related materials with magnetic ground states, which is ideal 
for the angle-resolved photoemission spectroscopy (ARPES) 
measurements. Subsequently, several groups have reported 
the direct observation of the bulk Weyl points and the surface 
Fermi arcs by ARPES in TaAs ll26H29ll , NbAs |[30l and TaP 
13111321 . confirming these compounds as TWSMs. For NbP, 
although a negative magnetoresistance induced by the chiral 
anomaly had been observed |[33l, an unambiguous evidence 
about the topological feature is still lack. 

On the other hand, NbP exhibit an ultrahigh carrier mobility 
of 5 X 10^ cm^ V“^ s“^ at 2 K |[33l|34l which is comparable 
to that of Cd 3 As 2 (9 x 10^ cm^ s”! at 5 K) |[35l. Un¬ 
like other TDSMs and TWSMs whose carries are either elec¬ 
tron or hole, NbP naturally hosts both types of carriers. Con¬ 
sequently, excellent electron-hole compensation leads to the 
huge magnetoresistance of 850000% at 1.85 K llMl exceed¬ 
ing that of WTe 2 (452700% at 4.5 K) (361. These properties 
make NbP even more intriguing. 

Therefore, probing the SS of NbP is of significant impor¬ 
tance not only to verify the topological character of the band 
structure, but also to understand its novel surface transport. 
Here, we report the observation of the surface electronic struc¬ 
ture of NbP on (001) surface by ultraviolet ARPES. Combin¬ 
ing with theoretical calculations, we identify odd number of 
bands cross the Fermi level in a closed loop, proving the exis¬ 
tence of unclosed FS originated from nontrivial surface states. 
The detailed connectivity of the Fermi arc is given by the cal¬ 
culations, which is qualitatively consistent with our experi¬ 


ment. Furthermore, the SSs exhibit unusual circular dichro- 
ism (CD) pattern, which might partially reflected the intrinsic 
information of the spin texture. Our experimental data and 
the first principle calculations not only establish that NbP is a 
TWSM, but also lay foundations for further investigations. 

High quality NbP single crystals were synthesized by two- 
step vapor transport technique as described elsewhere ( 23 . 
The crystal structure consists of alternate Nb and P layers, in 
which the adjacent NbP layers are rotated by 90° and shifted 
by a/ 2 , leading to a lack of inversion symmetry which is cru¬ 
cial to the realization of the TWSM [Fig. [Ja)]. Fig. [^b) 
presents the X-ray diffraction (XRD) pattern of our single 
crystal. The sharp diffraction peaks guarantee the high crys¬ 
talline quality of the samples. The typical dimension of 
the single crystal is about 1mm x 1mm x 0.5mm [inset in 
Fig.[2c)]. Resistivity data [Fig.[2c)] show that NbP is metal¬ 
lic over a wide temperature range (2 K to 300 K), while the 
residual resistivity ratio [RRR = p(300 K)/p(2 K)] is typically 
about 100. The predicted twelve pairs of Weyl points are il¬ 
lustrated in Fig. [Jd). For convenience, we denote the Weyl 
points in the = 0 plane as Weyl point 1 (Wl), and the others 
are denoted as Weyl point 2 (W2). 

ARPES measurements were performed at Beamline 105- 
ARPES of Diamond Light Source in which VG-Scienta 
R4000 electron analyzer is equipped. The typical angular res¬ 
olution is 0 . 2 ° and the overall energy resolution is better than 
20 meV. Samples were cleaved in-situ between adjacent NbP 
layers along the ( 001 ) plane [Fig.[ 2 a)], leaving fresh surface 
which violates C 4 rotation symmetry. The measurement pres¬ 
sure was kept below 9x10“^^ torr. 
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FIG. 3: (Color online) Observation of the Fermi arc SSs. (a) Schematic of the method to confirm the existence of the Fermi arc. Assuming that 
the FSs consist of one Fermi arc and one closed pocket (assigned as FSl and FS2, respectively). Then, choose a closed loop enclosing an odd 
number of Weyl points (green loop). If there exist a Fermi arc, an odd number of the Fermi crossings in the reference loop (3 times in this case) 
should be expected. The red and blue points represent the Weyl point with opposite chirality, (b) The illustration of the reference loop M - X 
- r - M with 3 Weyl points enclosed because the projection of W2 is two-fold degenerate. Note, the separation of Weyl points is exaggerated 
for clarity, (c) The measured band dispersions along high symmetry directions M - X - F. The upper panel shows the momentum distribution 
curves (MDCs) integrated near the Fermi level along corresponding high symmetry directions. The peak indicated by the red arrow in the 
MDCs along M - X - M corresponds to the residual spectra weight of the bands beneath Ep. (d) The calculated band dispersions along high 
symmetry directions M - X - F, showing excellent correspondence with the experimental data, (e) Calculated FS geometry around the X. The 
detailed evolution of different FS segments around the Weyl point is magnified in the inset. The size of the points in the calculated FS and 
band dispersions scale with the weight projected on the first unit cell of the P-terminated. (f) The measured FS showing overall agreement 
with the calculations, (g) The measured dispersions along cut 1 to 4 illustrated in panel (f). The corresponding EDCs are integrated around kp 
indicated by arrow in each cut. Solid lines are guide of eye for band dispersions. 


Our first principle calculation is performed by using the Vi¬ 
enna ah initio simulation package code 177113^ . An energy 
cutoff of 400 eV was adopted for the plane-wave expansion 
of the electronic wave function and the energy convergence 
criteria was set to 10“^ eV. Appropriate k-point meshes of (11 
X 11 X 1) and (100 x 100 x 1) were used for self-consistent 
and Fermi surface (FS) calculations, respectively. A k-point 
meshes of (400x400x1) were used to calculate the fine FS 
around Weyl points. The results are obtained by using the 
generalized gradient approximation Perdew-Becke-Erzenhof 
function |[^ . The spin-orbit coupling is taken into account 
by the second variation method 1401 . To simulate the (001) 
surface of NbP, we build a slab model composed of a seven 
unit cell, in which the top and bottom surface are terminated 
by P and Nb, respectively. A vacuum spacing of 20 A is used 
so that the interaction in the non-periodic directions can be 
neglected. The experimental lattice parameters were used and 


the positions of all the atoms were fully optimized. We pro¬ 
jected the band structure and FS to the first unit cell of the 
P-terminated side, in which the size of the red points scales 
with the projection weight. The obtained results fit the exper¬ 
imental data very well. 

The overall photoemission data of NbP are presented in 
Fig.[^ The measured FS consists of two parts centered around 
X and F, respectively [Fig.j^a)]. In general, each part is com¬ 
posed of two sets of orthogonal paddle-like FSs with some 
fine features inside, which is nearly in full accordance with 
the calculations [Fig.j^b)]. Figs.[^c) and (d) show the three- 
dimensional band structure around X and F, respectively. The 
measured band dispersions along F - M - A are again in ex¬ 
cellent agreement with the calculations. Here we note that not 
only the FSs are not four-fold symmetric, but also the band 
dispersions are anisotropic along X - M and Y - M. Since 
the cleaved (001) surface breaks the C4 rotational symme- 
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try, these results suggest that the observed band structure may 
originated from the SSs. The surface origin of the bands can 
be further confirmed by the photon-energy-dependent ARPES 
measurements in Fig. [^g), in which no evident dispersion 
along kz can be seen. The intensity variation of some peaks 
may be due to the matrix element effect. In fact, our pho¬ 
ton energy dependence results show that the band structure 
is dominated by these SSs even in a wider range (90 eV 
- 220 eV) (not shown here). This may be understood by the 
intrinsic surface sensitivity of ARPES due to the short escape 
depth of photoelectrons excited by vacuum ultraviolet pho¬ 
tons. 

One key signature of a TWSM is the existence of Fermi 
arcs originated from the topological SSs (H. To verify its ex¬ 
istence, we can choose a closed loop in the surface Brillouin 
zone with odd number of Weyl points enclosed and count the 
total time of Fermi crossings. As illustrated in Fig. [^a), the 
FS can only cross the loop an even number of times unless 
it is unclosed (Fermi arc). Although theoretically feasible, it 
is challenging to experimentally obtain the correct number of 
crossings in TaAs class of materials due to the extreme prox¬ 
imity of the W1 points and the degeneracy of the bands . 
To settle this problem, experimental data are commonly com¬ 
bined with calculations to help to precisely identify the Fermi 
crossings as exemplified in ref (7711 and (291 . In the present 
case, we choose M-X-F-Msls reference loop [Fig. [^b)] 
with 3 Weyl points included. The measured band dispersions 
are presented in Fig. |^c), which is in excellent correspon¬ 
dence with the calculations [Fig. [^d)]. We can identify 2, 5 
and 0 crossings with the assistance of our calculations along 
M - X, X - r and F - M, respectively. Thus, the surface FSs 
cross the reference loop 7 times, providing direct evidence of 
the existence of Fermi arcs on the (001) surface of NbP. 

Compared with the separation of the W1 points, that of the 
W2 points are larger. To examine the detail connectivity of 
the Fermi arc, we zoom in around the W2 Weyl point to check 
the SSs in higher resolution [Figs.|^(e) - (g)]. In Fig.j^e), 
our calculations show that the paddle-like FSs are actually 
composed of two dominant oval FSs. The inner one forms 
closed FS contributed by trivial SS. The outer one consists of 
three segments (Fermi arcs), two of which connect the W1 
and W2 points, while the rest one directly connects the W2 
points [Fig.[^(b)]. The Fermi arc connectivity here resembles 
that observed in NbAs (30ll . Fimited by the experimental res¬ 
olution, the measured FS only exhibits one broad oval feature 
[FigJ^f)]. To prove the correspondence with calculations, in 
Fig.j^g), we check the band evolution around the neck of the 
paddle-like FSs. From cut 1 to 4 [Fig.j^g)], a barely-resolved 
band (in red) can gradually be separated from the outer one 
(in yellow). This can be further verified by the energy distri¬ 
bution curves (EDCs) integrated around kf in each cut. The 
broad EDCs peak in cut 1 and 2 spilt into 2 peaks in cut 3 and 
4, suggesting that the broad feature is indeed contributed by 
two bands. The inner band further sinks below Ef forming 
closed FS, which is consistent with our calculations. Mean¬ 
while, the outer band further extends toward the W1 points. 



FIG. 4: (Color online) Circular dichroism of the electronic structure, 
(a) The Differential map of the RCP and LCP photoemission intensi¬ 
ties taken with 90 eV photons, (b) The corresponding intensity plot 
along several cuts indicated in panel (a). 


Although an open Fermi arc is not observed experimentally 
as them connects W1 and W2, the excellent correspondence 
between experiment and calculations establishes the existence 
of Fermi arc SSs in NbP. 

In TaAs class of materials which lack inversion symmetry, 
calculations show that Rashba-type spin polarization naturally 
exists in the the SSs and leads to strong spin textures (29l . 
With the presence of spin-orbit coupling, both spin and or¬ 
bital angular momentum (0AM) of a state would exhibit con¬ 
jugate textures around the FS. To further explore the spin tex¬ 
ture of the SSs, we performed CD-APRES experiment which 
has been demonstrated to be a powerful tool to probe the 0AM 
texture of the surface states in TIs (41114211 . 

Fig. I^a) show the FS mapping of the difference under 
the right-circular polarized (RCP) and left-circular polarized 
(TCP) light. The intensity map presents a strong difference 
between the RCP and the LCP data, which can be further ver¬ 
ified by several photoemission intensity plots along cut 5 to 
8. Around the Y point, the intensity of the oval FS sudden 
change its sign. This suggests that the CD observed here is 
most likely due to different 0AM, rather than some experi¬ 
mental geometry effects. Overall, the CD pattern is antisym¬ 
metric with respect to F, which can be explained by the time 
reversal symmetry. Due to the enormous amount of calcula- 
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tion needed, we leave the quantitative analysis of the CD to 
the future. 

To summarize, we have observed SSs on the (001) surface 
of NbP by using ultraviolet ARPES. Combined with calcu¬ 
lations, we identify the existence of the Fermi arc originated 
from topological SSs. Furthermore, we present the CD of the 
SSs, which may reflect the intrinsic information about the spin 
texture. Our results not only provide critical evidence to prove 
that NbP is a TWSM, but also lay foundation for further in¬ 
vestigations. 
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